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Abstract. A class of causal variational principles on a compact manifold is intro- 
duced and analyzed both numerically and analytically. It is proved under general 
assumptions that the support of a minimizing measure is either completely timelike, 
or it is singular in the sense that its interior is empty. In the examples of the circle, 
the sphere and certain flag manifolds, the general results are supplemented by a 
more detailed and explicit analysis of the minimizers. On the sphere, the minimal 
action is estimated from above and below. 
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1. Introduction 

Causal variational principles were proposed in [4J as an approach for formulating 
relativistic quantum field theory (for surveys see [71 [9]). More recently, they were intro- 
duced in a broader mathematical context as a class of nonlinear variational principles 
defined on measure spaces [8]. Except for the examples and general existence results 
in [6l [3l [8] and the symmetry breaking effect in the discrete setting [5j , almost nothing 
is known on the structure of the minimizers. In the present paper, we turn atten- 
tion to the analysis of minimizing measures in the continuous setting. For simplicity, 
we restrict attention to variational principles on a compact manifold which generalize 
the causal variational principles in spin dimension one with two prescribed eigenvalues 
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(see [SI Chapter 1] ) . But our methods are developed with a view to possible extensions 
to the non-compact setting and to a general spin dimension. 

More precisely, in Section [2] we introduce a class of causal variational principles on a 
compact manifold 3" and explain how this setting fits into the general context. As more 
specific model examples, we introduce variational principles on the circle, on the sphere, 
and on the flag manifold 3"^'^(C-^). In Section [3l we present numerical results on the 
sphere (see Figure [3]) and discuss all the main effects which will be treated analytically 
later on. In Section U we derive general results on the structure of the minimizers. We 
first derive the corresponding Euler-Lagrange equations and conditions for minimality 
(see Lemma 14.41 and Lemma 14. 5p . We then prove under general assumptions that the 
minimizers are either generically timelike (see Definition 14. 7p or the support of the 
minimizing measure p defined by 



is singular in the sense that its interior is empty (see Theorem 14.151 and Theorem 14. 17p . 
In the following sections, we apply these general results to our model examples and 
derive more detailed information on the minimizers. In Section [5l we consider the vari- 
ational principle on the circle. After briefly discussing numerical results (see Figured]), 
we prove a "phase transition" between generically timelike minimizers and minimizers 
with singular support and construct many minimizers in closed form (see Corollary 15. 21 
and Theorem 15. 4p . In Section [U the variational principle on the sphere is considered. 
We again prove the above phase transition (see Corollary 16.1 p and estimate the action 
from above and below (see Figure [6] and Proposition 16. 3p . Finally, in Section [7] we 
apply our general results to flag manifolds (see Theorem 17. ip . Moreover, we prove 
that minimizers with singular support exist (see Theorem 17. 2p and give an outlook on 
generically timelike minimizers. 



Before introducing our mathematical framework, we briefly put it in the general con- 
text. Causal variational principles can be formulated either in indefinite inner product 
spaces on an underlying space-time (the "space-time representation" ; see [U |6] and [El 
Chapters 3 and 4]) or in terms of the so-called local correlation matrices acting on 
the space of occupied particle states (the "particle representation" ; see [8l Chapters 1 
and 2]). Here we shall always work in the particle representation, whereas for the 
connection to the space-time representation we refer to the constructions in ^ Sec- 
tion 3.2] and [10]. Thus as in [Sj Chapters 1 and 2], we begin with a positive measure 
space (M, //) normalized to p{M) = 1. Moreover, for a given integer parameter /, 
we consider a closed subset 3" of the Hermitian (/ x /)-matrices. For technical sim- 
plicity, we here assume that 3" is a compact suhmanifold of Mat(C-^); this covers the 
variational principle with prescribed eigenvalues as considered in [8l Chapter 1 and 
Section 2.1]. Generally speaking, in a causal variational principle one minimizes a 
certain action S[F] under variations of a measurable function F : M — 3", imposing 
suitable constraints. Introducing the measure /9 on 3" by p{^) = /i(-F~^(r2)) (in other 
words, p = F^jj, is the push-forward measure), the action can be expressed by integrals 
over 3" (see [H Section 1.2]), 



supp p = {x £ 3' \ p{U) 7^ for every open neighborhood U of x} 



2. Preliminaries 




(2.1) 
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where the Lagrangian £ € C'^'^(3" x 9", M(j') is a given function. Now the only signifi- 
cance of the measure space (M, fi) is that it poses conditions on the possible form of 
the measure p. For example, in the discrete setting one chooses fi as the normalized 
counting measure on M = {1, . . . , m}; then the support of p necessarily consists of at 
most m points. However, in the continuous setting under consideration here, we do not 
want to impose any conditions on the measure p, but instead p should be allowed to be 
any normalized positive regular Borel measure on 3". Then the measure space (M, p,) 
is no longer needed. For simplicity, we also leave out additional constraints (like the 
trace or identity constraints; see [SJ Section 2.3]). This leads us to the following setting: 
Let 3" be a smooth compact manifold (of arbitrary dimension). For a given function 



(2.2) 
(2.3) 

(2.4) 

(2.5) 



X> G C°°(J X 3", M) being symmetric: D(x, y) = Ti{y, x) V x, y G J 

and strictly positive on the diagonal: T>(x,x) > 

we define the Lagrangian C by 

£ = max(0, V) G C°'^( J x J, M+) . 

Introducing the action S by (|2.ip . our action principle is to 

minimize S under variations of p G , 

where 9Jt denotes the set of all normalized positive regular Borel measures on 3". In 
view of the symmetric form of (|2.ip . it is no loss of generality to assume that C{x,y) is 
symmetric in x and y. Therefore, it is natural to assume that also D{x,y) is symmet- 
ric (j2.2p . If (j2.3p were violated, every measure supported in the set {x : T){x,x) < 0} 
would be a minimizer. Thus the condition (j2.3p rules out trivial cases. 

The existence of minimizers follows immediately from abstract compactness argu- 
ments (see [8, Section 1.2]). 

Theorem 2.1. The infimum of the variational principle (j2.5p is attained in 9Jt. 

Note that the minimizers will in general not be unique. Moreover, the abstract frame- 
work gives no information on how the minimizers look like. 

The notion of causality can now be introduced via the sign of 2). 

Definition 2.2 (causal structure). 

timelike 
lightlike 
spacelike 



Two points x,y £ 3' are called < 



separated if < 



D{x,y) > 
D(x,2/) = 
D{x,y) < . 



We define the sets 



= {y G 3" with T){x,y) > 0} 
= {y G 3" with D(x,y) > 0} 
= dl{x)nd{3'\j{x)) 

if and only if the function D(x, .) 



open lightcone 
closed lightcone 
boundary of the lightcone . 

changes sign in every neighborhood 



I{x) 
J{x) 
lC{x) 

Thus y G /C(x) 
of y. 

Our action is compatible with the causal structure in the sense that if x and y have 
lightlike or spacelike separation, then the Lagrangian vanishes, so that the pair [x, y) 
does not contribute to the action. Note that for a given minimizer p, we have similarly 
a causal structure on its support. 

In order to work in more specific examples, we shall consider the following three 
model problems. 
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Figure 1. The function D. 



(a) Variational principles on the sphere: 

We consider the setting of [H Chapter 1] in the case / = 2 (see also [8^ Exam- 
ples 1.5, 1.6 and 2.8]). Thus for a given parameter r > 1, we let 3" be the space 
of Hermitian (2 x 2)-matrices whose eigenvalues are equal to 1 + r and 1 — r. 
Writing a matrix F G 3" as a linear combination of Pauli matrices, 

F = T x-a + 1 with x G S"^ C , 

we can describe F by the unit vector x (here • denotes the scalar product 
in R^). Thus 3" can be identified with the unit sphere S"^. The function D is 
computed in [HI Example 2.8] to be 

D(x, y) = 2r2 (1 + (x, y)) (2 - (1 - (x, y))) . (2.6) 

This function depends only on the angle "dxy between the points x,y ^ S"^ 
defined by cos-dxy = {x,y)- Considered as function of •& G [0, vr], T> has its 
maximum at = and is minimal if cos('!?) = — r^^. Moreover, T>{tt) = 0. 
Typical plots are shown in Figure [H In the case r > 1, the function D has two 
zeros at vr and 

■i^max := arccos ^1 - . (2.7) 

In view of ()2.4p . the Lagrangian is positive if and only if < "i? < ??max- 
Thus X(x) is an open spherical cap, and J{x) is its closure together with the 
antipodal point of x, 

X(x) = {y: (x,y) >1-^}, J{x) =T{^ U {-x} 

If T increases, the opening angle t?max of the lightcones gets smaller. In the 
degenerate case r = 1, the function D is decreasing, non-negative and has 
exactly one zero at = vr. Hence the Lagrangian C coincides with D. All 
points on the sphere are timelike separated except for antipodal points. The 
lightcones are T(x) = S'^\{—x} and J7'(x) = 5^. 

If we regard p as a density on the sphere, the action (12. ip looks like the energy 
functional corresponding to a pair potential C (see for example |I3]). Using 
physical notions, our pair potential is repelling (because C{i)) is a decreasing 
function) and has short range (because C vanishes ii i} > i?max)- 
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(b) Variational principles on the circle: 

In order to simplify the previous example, we set 3" = S^. For D we again 
choose (j2.6p . 

(c) Variational principles on the flag manifold 3^^''^{C^): 

As in [8], Chapter 1] , for a given parameter r > 1 and integer parameters / > 2 
we let 3" be the space of Hermitian (/ x /)-matrices of rank two, whose nontrivial 
eigenvalues are equal to 1 + r and 1 — r. Every x G 3" is uniquely described by 
the corresponding eigenspaces U and V. By considering the chain U C (UUV), 
X can be identified with an element of the flag manifold 3^^'^(C-^), the space 
of one-dimensional subspaces contained in a two-dimensional subspace of C-^ 
(see It is a (4/ — 6)-dimensional compact manifold. Every U G U(/) gives 

rise to the mapping x — t- UxU~^ on 3". This resulting group action of U(/) 
on 3" acts transitively, making 3" to a homogeneous space (see [11] for details) . 

For two points x,y G 3", we denote the two non-trivial eigenvalues of the 
matrix product xy by , \^ G C and define the Lagrangian by 

c{x,y) = \{\xiy\-mf . 

This Lagrangian is U(/)-invariant. In order to bring it into a more convenient 
form, we first note that by restricting to the image of y, the characteristic 
polynomial of xy changes only by irrelevant factors of A, 

det(xy - Al) = A-^~^ det {{-Kyxy - Al)|imj/) , 

where iTy denotes the orthogonal projection to Imy. It follows that \^ and \^ 
are the eigenvalues of the (2 x 2)-matrix ■Kyxy\ijay In particular, 

X'^yX'^y = d(ii{'KyXy\lray) = d(ii{-KyXT:y\lray) det(y|imj/) > , 

because the operator TTyXiTy again has at most one positive and one negative 
eigenvalue. Moreover, the relation X^ + Xj^ = Tr(xy) G M shows that the 
two eigenvalues are either both real and have the same sign or else form a 
complex conjugate pair, in which case the Lagrangian vanishes. Finally, us- 
ing that (A+^) + {Xj') = Tr((xi/)2), the La grangian can be written in the 
form with 

D(x,y) = i (A:;^ - A^^)' = Tr{{xyf) - ^(TY(xy))' . (2.8) 

We finally comment on the limitations of our setting and mention possible general- 
izations. First, we point out that our structural results do not immediately apply 
in the cases when 3" is non-compact or when additional constraints are considered 
(see [H Chapter 2]). However, it seems that in the non-compact case, our methods 
and results could be adapted to the so-called moment measures as introduced in [51 
Section 2.3]. A promising strategy to handle additional constraints would be to first 
derive the corresponding Euler-Lagrange equations, treating the constraints with La- 
grange multipliers. Then one could try to recover these Euler-Lagrange equations as 
those corresponding to an unconstrained variational problem on a submanifold S C 3", 
where our methods could again be used. We finally point out that in the case of 
higher spin dimension n > 1, it is in general impossible to write the Lagrangian in 
the form (j2.4p with a smooth function D, because the Lagrangian is in general only 
Lipschitz continuous in the open light cone. A possible strategy would be to first 
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show that the support of p hes on a submanifold S C 3", and then to verify that by 
restricting £ to 9 x S, it becomes smooth in the open hghtcones. 

3. Motivation: Numerical Results on the Sphere 

In order to motivate our general structural results, we now describe our findings in 
a numerical analysis of the variational principle on the sphere (see Example (a) on 
paged]). Clearly, in a numerical study one must work with discrete configurations. 
Our first attempt is to choose a finite number of points xi, . . . ,Xm £ 5"^ and to let p 
be the corresponding normalized counting measure, i.e. 

m 

fdp:=-S2f{xi) yf€C^{S^). (3.1) 

i=l 

Then the action (|2.ip becomes 

^ m 

^=^Y.^(^^''^jy (3.2) 
III' • • , 

By varying the points Xj for fixed m, we obtain a minimizer pm- Since every normalized 
positive regular Borel measure can be approximated by such counting measures, we 
can expect that if we choose m sufficiently large, the measure pm should be a good 
approximation of a minimizing measure p € dJl (more precisely, we even know that 
Pm ^ P as m — oo with convergence in the weak (C'^)*-topology). 

If r is sufficiently large, the opening angle of the lightcones is so small that the m 
points can be distributed on the sphere such that any two different points are spacelike 
separated. In this case, the action becomes 

5 = —£(?? = 0) , 

m 

and in view of ()3.2p this is indeed minimal. The question for which r such a configu- 
ration exists leads us to the Tammes problem, a packing problem where the points are 
distributed on the sphere such that the minimal distance between distinct points is 
maximized, see jl4| . More precisely, we know that the Tammes distribution is a mini- 
mizer of our action if r is so large that 'dm > "f^max- Until now, the Tammes problem 
is only solved if m < 12 and for m = 24 (for details see [2j and the references therein). 
For special values of m, the solutions of the Tammes problem are symmetric solids 
like the tetrahedron (m = 4), the octahedron (m = 6), the icosahedron (m = 12) and 
the snub cube (m = 24). Moreover, much research has been done on the numerical 
evaluation of spherical codes, mostly by N.J. A. Sloane, with the collaboration of R.H. 
Hardin, W.D. Smith and others, |14) . containing numerical solutions of the Tammes 
problem for up to 130 points. 

In the case "i^m < 'dmax, the measure pm was constructed numerically using a sim- 
ulated annealing algorithn£|. In order to get optimal results, we used this algorithm 
iteratively, using either a Tammes distribution or previous numerical distributions 
as starting values. Using that 2) depends smoothly on r, it is useful to increase or 
decrease r in small steps, and to use the numerical minimizer as the starting con- 
figuration of the next step. In Figure [21 the numerically found is plotted for 
different values of m as a function of the parameter r. The resulting plots look rather 

"'^We use the "general simulated annealing algorithm" by J. Vandekerckhove, © 2006, 
http : //www.mathworks .de/matlabcentral/fileexchange/ 10548. 
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Figure 2. Numerical minima for the counting measure on the sphere. 

compHcated. The considered values for m are too small for extrapolating the limiting 
behavior as m — )■ oo. Nevertheless, one observation turned out to be very helpful: 
Near r ~ 1.2, the plots for different values of m look the same. The reason is that 
some of the Xi coincide, forming "clusters" of several points. For example, in the 
case m = 12, the support of p only consists of six distinct points, each occupied by 
two Xi. A similar "clustering effect" also occurs for higher r if m is sufficiently large. 

These findings give the hope that for large m, the minimizers might be well- 
approximated by a measure supported at a few cluster points, with weights count- 
ing the number of points at each cluster. This was our motivation for considering a 
weighted counting measure. Thus for any fixed m, we choose points xi, . . . ,Xm S 5*^ 
and corresponding weights pi, . . . , pm with 

m 

Pi > and Pi = 1 . 

1=1 

We introduce the corresponding measure p in generalization of (j3.ip by 

„ m 

/ fdp ■.= Y^p^f{xi) V/gC°(52). (3.3) 
•^^^ i=i 

Seeking for numerical minimizers by varying both the points Xi and the weights pi, we 
obtain the plots shown in Figure El 

These plots suggest the following structure of the minimizers. Let us denote the the 
minimizing weighted counting measure for a given m by p{m). Then for any fixed r, 
the series S[p{m)\ is monotone decreasing (this is obvious because every p{m) can be 
realized by a weighted counting measure with rrij^ > m summands by choosing m+ — m 
weights equal to zero). The important observation is that there is an integer mo from 
where on the series stays constant, i.e. 

S[p{m^)] > S[p{mo)] = 5[p(m+)] V m_ < niQ < m+ . 

This implies that the measure is also a minimizer in the class of all Borel measures. 
This leads us to the following 



8 



F. FINSTER AND D. SCHIEFENEDER 




X 



Figure 3. Numerical minima for the weighted counting measure on the sphere. 

Conjecture 3.1. For any r > 1, there is a minimizer p £ dJl of the variational 
problem on the sphere which is a weighted counting measure supported at mo points. 

Prom Figure [3] we can read off the value of mo as a function of r. Generally speaking, 
mo increases as r gets larger. This corresponds to the fact that for increasing r, the 
opening angle i?max of the light cones gets smaller, so that it becomes possible to 
distribute more points on the sphere which are all spatially separated from all the 
other points. 

The more detailed numerical study of the minimizers showed another interesting 
effect. For values r < Tc := \/2, we found many different minimizers of different form. 
They all have the property that they are completely timelike in the sense that all points 
in the support of the minimizing measure have timelike or lightlike separation from all 
the other points. We found minimizers supported on an arbitrarily large number of 
points. If on the other hand r > Tc, all minimizers were supported on at most mo(T) 
points, indicating that every minimizing measure p G 9Jt should be discrete with finite 
support. The intermediate value t = Tc correspond to the opening angle i^max = f of 
the light cones. 

Conjecture 3.2. If t < Tc, every minimizer is completely timelike. If conversely 
T > Tc, every minimizing measure is discrete with finite support. 

More graphically, one can say that for r > Tc, our variational principle spontaneously 
generates a discrete structure on the sphere. The two regions r < Tc and t > can 
also be understood as two different phases of the system, so that at r = Tc we have a 
phase transition from the completely timelike phase to the discrete phase. 

The above numerical results will be the guide line for our analysis. More precisely, 
the completely timelike phase will be analyzed in Section 14.21 using the notion of 
"generically timelike", whereas in Section S]3] we will develop under which assumptions 
and in which sense the support of the minimizing measure is discrete or "singular" . The 
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phase transition is made precise in Theorem 14.151 and 14.171 by stating that minimizing 
measures are either generically timehke or singular. 

4. General Structural Results 

We now return to the general variational principle (j2.5p with the Lagrangian of the 
form (|2.4p and (|2.2p on a general smooth compact manifold 3". Let us introduce some 
notation. For a given measure /> G 9Jt we define the functions 



£ix) = j^C{x,y)dp{y) ^ C''^\y) (4.1) 
d(x)= / V{x,y)dp{y) eC~(J). (4.2) 



Moreover, we denote the Hilbert space L?'{'J,dp) by (;Kp,(.,.)p) and introduce the 
operators 



Cp : ^p^'Kp : {Cp^){x) = / C{x,y) ij{y) dp{y) (4.3) 

Dp : 'Kp^'Kp -.^^ {'J)pij){x) = [ V{x,y) ^{y) dp{y) . (4.4) 



Lemma 4.1. The operators Cp and Dp are self-adjoint and Hilbert- Schmidt. The 
eigenf unctions of Cp ( and Dp) corresponding to the non-zero eigenvalues can he ex- 
tended to Lipschitz-continuous (respectively smooth) functions on 3". 

Proof. We only consider Dp, as the proof for Cp is analogous. The self-adjointness 
follows immediately from the fact that T){x,y) is symmetric. Moreover, as the kernel 
is smooth and 3~ is compact, we know that 

\'D{x,y)\'^dp{x)dp{y) < oo . 

This implies that Dp is Hilbert-Schmidt (see \12\ Theorem 2 in Section 16.1]). 
Suppose that Dptp = Xtp with A / 0. Then the representation 

i^i^) = T I '^ix,y)il){y)dp{y) 



shows that ^ G C°°(3'). □ 

The following notions characterize properties of 3~ and the function D which will be 
needed later on. 

Definition 4.2. A measure p & dJl is a homogenenizer o/D ifsuppp = 3" and both 
functions 



ifiix) := / C{x,y)dp{y) and d^(x) := / D{x,y) dp{y) 

are constant on 3". The function D is called homogenizable if a homogenizer exists. 

In the Examples (a), (b) in Section [21 we can always choose the standard normalized 
volume measure as the homogenizer. More generally, in Example (c) we choose for p 
the normalized Haar measure, obtained by introducing a U(/)-invariant metric on 3" 
and taking the corresponding volume form (see for example [l] Section L5]). 

The next proposition gives a sufficient condition for a homogenizer to be a minimizer. 
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Proposition 4.3. // > 0, the homogenizer fi is a minimizer of the variational 
principle (j2.5p . 

Proof. We denote the constant function on 3" by I3- = 1. If is a homogenizer, 
this function is an eigenfunction of which can be completed to an orthonormal 
eigenvector basis (■i/'OieNo of '^ith "00 = Ij and corresponding eigenvalues Aj > 0. 
Using an approximation argument in the C'^(J-")*-topology, it suffices to show that 

for any tp G (7*^(3") with "0^0 and {ip, Ig-)^ = 1. We write ip in the eigenvector 
basis 0i, 

00 

ip = y^^citpi . 

The condition {tp, Is^)^ = 1 implies that cq = 1. Thus 

00 

S[ipij] = {ip,C^ip)^ = Ao + ^ IcjpAi > Ao = S[fj] . |-| 

1=1 

4.1. The Euler-Lagrange Equations. Let us assume that p is a minimizer of the 
variational principle (|2.5|) . 

S[p] = inf S[p]=: Smin ■ 

We now derive consequences of the minimality. In the first lemma, we consider first 
variations of p to obtain the Euler-Lagrange equations corresponding to our variational 
principle. The second lemma, on the other hand, accounts for a nonlinear effect. 

Lemma 4.4. (The Euler-Lagrange equations) 

^Isuppp = ^ — '^min • 

Proof. Comparing ()2.1|) with ()4.ip . one sees that 

5min = [ idp. (4.5) 

Since i is continuous and 3" is compact, there clearly is y G 3" with 

e{y) = mie. (4.6) 

We consider for t G [0, 1] the family of measures 

Pt = {l-t)p + t5y em, 

where 6y denotes the Dirac measure at y. Substituting this formula in (j2.ip and 
differentiating, we obtain for the first variation the formula 

6S := lim = _25^,, + 2£(y) . 

Since p is a minimizer, 6S is non-negative. Combining this result with ()4.5p and ()4.6p . 
we obtain the relations 

mie = £{y) > Smin= [idp. 
It follows that i is constant on the support of p, giving the result. □ 
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Lemma 4.5. The operator Cp is non-negative. 
Proof. Lemma (j4.4p yields that for any x £ supp p, 

(£plj)(x) = j^C{x,y) dp{y) = = 5min ^7{x) , 

showing that the constant function Ig- is an eigenvector corresponding to the eigen- 
value 5min > 0. 

Assume that the lemma is wrong. Then, as Cp is a compact and self-adjoint operator 
(see Lemma [4.ip . there exists an eigenvector ijj corresponding to a negative eigenvalue, 
Cpip = Xip with A < 0. We consider the family of measures 

In view of Lemma 14.11 is continuous and therefore bounded. Thus for sufficiently 
small \t\, the measure pt is positive. Moreover, the orthogonality of the eigenfunc- 
tions Ij and ip implies that 

Pt{J) = J^l^{lj + tiP)dp=l + t{l:f,iP)p = l, 

showing that pt is again normalized. Finally, again using the orthogonality, 

5[pt] = {{l^ + tlP),Lp{lj + tlP))p = Smin + Xt^ (V','0>p • 

Thus Pt is an admissible variation which decreases the action, a contradiction. □ 

An immediate consequence of this lemma is a useful positivity property of the La- 
grangian when evaluated on a finite number of points in the support of p. 

Corollary 4.6. For a finite family xq, . . . , X]\r £ supp/? (with N G Nj, the Gram 
matrix L defined by 

L = (C{xi,xj)) 
is symmetric and positive semi-definite. 

Proof. Given e > and a vector u = [uq, . . . ,un) G C^^^, we set 

N 

'f^^('^^ = E^(il^^XB4.d^)e:Kp, 

where is a ball of radius e (in a given coordinate system) . Lemma 14.51 implies 
that {ijje, Cpil)^) > 0. Taking the limit e \ 0, it follows that 

{u,Lu)cN+i = lim{lpe,Cp'4'e)p > . ^ 
e\0 U 

4.2. Generically Timelike Minimizers. 

Definition 4.7. A minimizing measure p G 9K is called generically timelike if the 

following conditions hold: 

(i) D{x,y) > for all x,y £ supp p. 

(ii) The function d defined by (j4.2p is constant on 3'. 

This constant can easily be computed. 

Lemma 4.8. Suppose that p is a generically timelike minimizer. Then 

d (x) = 5min for all X £ 3' . 
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Proof. Since £ and D coincide on the support of p, we know that 



5min = / / ^x, y) dp{x) dp{y) =11 y) dp{x) dp{y) . 

J J J J 7x7 

Carrying out one integral using ()4.2p . we obtain 

5min = / cl(x) dp{x) , 

J7 

giving the result. □ 

In the remainder of this subsection, we assume that D is homogenizable (see Defi- 
nition [32]) and denote the homogenizer hy p £ 9Jt. 



Lemma 4.9. IfD^ has only a finite number of negative eigenvalues, the kernel D{x,y) 
has the representation 

N 

Vix, y) = + ^ Cl)n{x) (t)n{y) (4.7) 
n=l 

with G N U {oo}, f„ G M, m„ 7^ 0, and 4>n £ C°°(3"), where in the case N = 00 the 
series converges uniformly. 

Proof. By definition of the homogenizer, the function Ig^ = 1 is an eigenfunction of 
the operator D^. Denoting the corresponding eigenvalue by uq, we obtain the spectral 
representation (|4.7p . 

If is positive semi-definite, the uniform convergence is an immediate general- 
ization of Mercer's theorem (see [121 Theorem 11 in Chapter 30], where we replace 
the interval [0, 1] by the compact space 3", and the Lebesgue measure by the mea- 
sure p). In the case when T>^ has a finite number of negative eigenvalues, we apply 
Mercer's theorem similarly to the operator with kernel T){x,y) — Yli^i '^i{x)'^i{y) ^ 
where \i,...Xk are the negative eigenvalues with corresponding eigenfunctions ipi. 
By construction, this operator is positive semi-definite, and in view of Lemma [4.1 1 its 
kernel is continuous. □ 

Lemma 4.10. Suppose that p is a generically timelike minimizer and that the opera- 
tor Dfj^ has only a finite number of negative eigenvalues. Then 

S[p] = uq and / (pniv) dp{y) = for all n e {I, . . . ,N] . 
J7 

Proof. Using the decomposition of the kernel (|4.7p and the uniform convergence, we 
obtain 

^ r 

d{x) = VQ + ^Un(t)n{x) I (t)n{y) dp{y) . 
n=l ''^ 

Applying Lemma [4.81 gives the claim. □ 
Proposition 4.11. Suppose that is a positive semi-definite operator on 'Ji^j^. Then 

In the case of equality, every minimizer is generically timelike. 
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Proof. If T>^ is positive semi-definite, all the parameters i^n in (|4.7p are positive. It 
follows that for every measure p G 9Jt, 

y) dp{x) dp{y) > [[ 1){x, y) dp{x) dp{y) > p{3^f = . (4.8) 
Jx5 JJsxS 



Let us assume that equality holds. It then follows from (j4.8p that C and D coincide 
on the support of p and thus 2)(x,y) > for all x,y £ suppp. Moreover, we find 
from (|^ that 

N 2 

= + X] / ^"^^y^ ' 

n=l 

and thus 

/ 4>n{y) dp = for all n > 1 . 

It follows that dp is a constant. We conclude that p is generically timelike. □ 

This proposition can be used to construct generically timelike minimizers. 

Corollary 4.12. Suppose that T>^ is a positive semi-definite operator on Jf^. Assume 
that the function f £ has the following properties: 

(a) Ti{x, y) = £(x, y) for all x,y £ supp /. 

(b) / f{x)dp{x) = l and I f{x)(l)n{x)dp{x)=0 for all n £ {1, . . . , N}. 
J? J J 

Then the measure dp = f dp is a generically timelike minimizer. 

Proof. The assumption (a) implies that 

S[p]= ff T>{x,y)dp{x)dp{y). 

Using the decomposition (|4.7|) and the relations (b), we find that S[p\ = vq. We now 
apply Proposition 14.111 □ 

We conclude this section by stating obstructions for the existence of generically 
timelike minimizers. 

Proposition 4.13. Assume that one of the following conditions hold: 

(I) The operator T>^ has only a finite number of negative eigenvalues, and the 

eigenvalue fQ in the decomposition (j4.7p is not positive. 
(II) For every x G 3" there is a point y £ 3' with J{x) n J{y) = ("condition of 
disjoint lightcones"). 

(Ill) For every x G 3" there is a point —x ^ T{x) with J{x) = I{x) U {— x} and 
I{x) nX(— x) = ( "condition of antipodal points"). 
Then there are no generically timelike minimizers. 

Proof. We first show that 5min > 0. Namely, choosing x in the support of a minimizing 
measure p, we know from ()2.3p and the continuity of T> that there is a neighborhood U 
of X and 5 > such that D{x, y) > 6 for all y £ U. It follows that 

C{x,y)dpix)dp{y)>5p{Uf>0. 

UxU 

Case (I) is obvious in view of Lemma 14.101 and the fact that 5min > 0. To prove 
the remaining cases (II) and (III) , we assume conversely that there exists a generically 
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timelike minimizer p G 9K. Choosing a point x G supp p, we know from property (i) in 
Definition 14. 71 that supp/o C J{x). In case (II), we choose y E 3" with J {x)r\J {y) = 
to obtain 



d(y)= / !D(y,z)dp(z) <0< 

JJ{x) 



in contradiction to Lemma |4.8[ 

In case (III), we know that suppp C J{x) = T{x) U {—x}. If —x ^ supp/?, the 
estimate 

d(-x) = [ V{-x,z)dp{z) = [ T»(-x,z)dp(z) < < 5mm 

Jj{x) Jx(x) 

again gives a contradiction, where in (*) we used that X{x) n X(— x) = 0. If con- 
versely —X G supp p, then supp p C J{x) n i7(— x) = {x} U {— x} (where we again used 
that X[x) nX(— x) = 0). Hence the integral in (j4.2p reduces to a sum over two points, 

d(y) = p{{x}) V{y, x) + p{{-x}) V{y, -x) . (4.9) 

In view of our assumption (j2.3p . we know that x G X(x). On the other hand, the 
relation X(x) nX(— x) = shows that — x ^ X(x). Hence there is a point y S dX{x). It 
follows that 2)(?/, x) = (because y G dl{x)) and also 2)(y, — x) < (because y £ I{x) 
and thus y ^ X(— x)). Using these inequalities in ()4.9p . we again find that d(y) < 0, a 
contradiction. □ 

It is an interesting question how the support of a generically timelike minimizer p may 
look like. The next proposition (which will not be used later on) quantifies that supp p 
must be sufficiently "spread out" . 

Proposition 4.14. Assume that p is a generically timelike minimizer and that the 
operator has only a finite number of negative eigenvalues. Then every real func- 
tion ijj G 2)^(5{^) with 



V'(x) dp{x) = (4.10) 
changes its sign on the support of p (here p is again the homogenizer of Definition\^.^. 



Proof. We return to the spectral decomposition (]4.7p of the operator 2)^. Since the 
eigenfunctions are orthogonal in Jf^, we know that 



bn dp = for all n > 1 . 
Representing ^jJ in an eigenvector basis of T>^ and using (|4.10p . we find 

N 



n=l 



with complex coefficients k„. Integrating with respect to p, we can apply Lemma 14.101 
to obtain 

TV 

/ '>p{x)dp{x) = y2Kn / 0„(x) dp(x) = . 
Hence ip changes its sign on the support of p. □ 
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4.3. Minimizers with Singular Support. We now state results on the support of 
a minimizing measure. 

Theorem 4.15. Let "J he a smooth compact manifold. Assume that T){x,y) is sym- 
metric (|2.2|) and equal to one on the diagonal, T>{x,x) = 1. Furthermore, we assume 
that for every x G 3" and y £ IC{x), there is a smooth curve c joining the points x 
and y, along which D{.,y) has a non-zero derivative at x, i.e. 

|D(c(t)..)|,^„#0, (4.11) 

where we parametrized the curve such that c(0) = x. Then the following statements 
are true: 

(A) // 3", T> are real analytic, then a minimizing measure p is either generically 
timelike or supp/o = 0. 

(B) IfD is smooth and if there is a differential operator A on C°°(3") which vanishes 
on the constant functions such that 

A^V{x, y) <0 for all y G I{x) , (4.12) 

o 

then supp p = 0. 

A typical example for A is the Laplacian corresponding to a Riemannian metric on 3". 
Note that the condition (|4.1ip implies that for every y G 3", the set {x | y G /C(x)} is a 
smooth hypersurface, which the curve c intersects transversely (in the applications of 
Section [5] and m this set will coincide with /C(y), but this does not need to be true in 
general) . 

The condition ()4.1ip can be removed if instead we make the following symmetry 
assumption. 

Definition 4.16. The function T) is called locally translation symmetric at x with 
respect to a curve c(t) with c(0) = x if there is e > and a function f G C°^((— 2e, 2e)) 
such that the curve c is defined on the interval (— e, e) and 

V{c{t),c{t')) = f{t - t') for all t, t' G {-e, s) . 

Theorem 4.17. Let 3' be a smooth compact manifold. Assume that T>{x,y) is sym- 
metric (j2.2p and strictly positive on the diagonal (j2.3p . Furthermore, we assume that 
for every x G 3~ and y G /C(x), there is a smooth curve c joining the points x and y 
such that T> is locally translation symmetric at x with respect to c, and such that the 
function T>{c{t),y) changes sign at t = (where we again parametrize the curve such 
that c(0) = x). Then statement (A) of Theorem \4-15\ holds, provided that the curve c 
is analytic in a neighborhood oft = 0. Assume furthermore that there is p £N with 

dP 

-V(c{t),y) /O. (4.13) 



dtp 



t=o 



Then statement (B) of Theorem \4 .151 again holds. 

In the smooth setting, the above theorems involve quite strong additional assump- 
tions (see (14. lip . (I4.12P and (14.131) ). The following counter example shows that some 
conditions of this type are necessary for the statements of these theorems to be tru^. 



9 

We would like to thank Robert Seiringer for pointing out a similar example to us. 
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Example 4.18. Let f,g G Cq°([— tt, vr]) be non-negative even functions with 

supp/c[-f,f] , SUPP5 C (-7r,-f] U [f,7r) . 
We introduce the function D G C°°(52 x S'^) by 

Vix, y) = -5(dist(x, y)) + [ /(dist(x, z)) /(dist(z, y)) dfi{z) , (4.14) 

where dfi is the standard volume measure, and dist denotes the geodesic distance 
(taking values in [0,7r]). Note that the two summands in (j4.14p have disjoint supports 
and thus the corresponding Lagrangian ()2.4p simply is 

Cix,y)= [ /(dist(x,z))/(dist(z,y))(i^(z), (4.15) 

We again consider T){x,y) and C{x,y) as the integral kernels of corresponding opera- 
tors T)^ and £^ on the Hilbert space IK^ = L?'{S'^, dfi). 

First, it is obvious that T)[x,y) is symmetric and constant on the diagonal. Next, 
it is clear by symmetry that the measure /i is a homogenizer (see Definition 14. 2p . 
Moreover, writing £^ as = where is the operator with integral kernel /, one 
sees that the operator is non-negative. Thus by Proposition 14.31 the measure is 
minimizing. If the function g is non-trivial, there are points x, y which are spacelike 
separated, so that this minimizer is not generically timelike. Also, its support obviously 
has a non- vanishing interior. We have thus found a minimizing measure which violates 
statement (A) of Theorem 14.151 

The remainder of this section is devoted to the proof of the above theorems. We 
begin with a simple but very useful consideration. Suppose that for given x G 9", the 
boundary of the light cone /C(x) does not intersect the support of p. As the support 
of p is compact, there is neighborhood {/ of x such that 

K.{z) n supp p = for all 2; G [/ . 

Thus introducing the measure p = Xx{x) the function £ can for all z £ U be repre- 
sented by 

i{z)= [ C[z,Odm= I D{z,Odm- (4.16) 

This identity can be used both in the smooth and in the analytic case. 

Lemma 4.19. // ()4.12p holds, then for every x G suppp the set /C(x) H supp/o is 
nonempty. 

Proof. Applying the Laplacian to ()4.16p gives 

AJ{x) = A^D(x, z) dp{z) < , 

where in the last step we used (j4.12p and the fact that x G supp p. This is a contra- 
diction to Lemma 14.41 □ 

Lemma 4.20. Suppose that 3" and D are real analytic. Assume that there exists 
a point X G suppp such that JC{x) Pi suppp = 0. Then p is generically timelike 
and supp p C X(x). 
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Proof. We introduce on 3" the function 

d(y) = / 'D{y,z)dp(y) 



Then d is real analytic and, according to (j4.16p . it coincides on U with the function (.. 
Since x G supp p, the Euler-Lagrange equations in Lemma (|4.4|) yield that ^ = 5min in 
a neighborhood of x. Hence d = 5min in a neighborhood of x, and the real analyticity 
implies that 

d = 5min on J . 

It follows that 

5min = / A{x)dp{x) = D{x,y) dp{x) dp{y) 

^^^^^ (4.17) 
< // ^x,y) dp{x) dp{y) = / £{x) dp{x) = S^in K^) , 

and thus p{F) = 1. Since p < p and p is normalized, we conclude that p = p. 
Thus d = d = 5min- Moreover, the inequality in (|4.17p becomes an equality, showing 
that £ = D on the support of p. Thus p is indeed generically timelike. □ 

To complete the proof of Theorems 14.151 and 14.171 it remains to show the following 
statement: 

K,{x) n supp/9 = for all x G suppp . (4-18) 

We proceed indirectly and assume that there is a point y G lC{x) n suppp. Our 
strategy is to choose points xq, . . . , in a neighborhood of x such that C restricted to 
the set {xo, . . . , x^, y} is not positive semi-definite, in contradiction to Corollary 14.61 
The points xq, . . . , x^ will all lie on a fixed smooth curve c which joins x and y chosen 
as in the statement of the theorems. We parametrize c such that c(0) = x and c(l) = y, 
and by extending the curve we can arrange (possibly by decreasing e) that the curve 
is defined on the interval {—ke, 1]. By the assumptions in Theorems 14.151 and 14.171 we 
know that T){c{t),y) changes sign at t = 0. Depending on the sign of D(c(e),0), we 
introduce the equidistant "chain" of points 

fxo = c(e), xi = c(0), X2 = c{-e), ... ,xk = c{-{k-l)e) if D(c(e), 0) > 

{ (4.19) 

[xo = c(— e), xi = c(0), X2 = c(e), ... , x^ = c((/c — l)e) if D(c(e), 0) < 0. 

(thus y has timelike separation from xq, lightlike separation from xi = x, and spacelike 
separation from X2, . . . ,Xfc). Then by construction, xq G ^{y), whereas all the other 
points of the chain are spacelike or lightlike separated from y. 

For the proof of Theorem I4.15| it suffices to consider a chain of three points. 

Lemma 4.21. Assume that T){x,y) is symmetric (j2.2p and equal to one on the diag- 
onal, D(x,x) = 1. Then for xq,xi,X2 as given by ()4.19p in the case k = 2, there is a 
real constant ai such that for all sufficiently small e, 

V{xi,Xj) = l + ai\i- £^ + 0{e^) for all i,j £ {0,1,2} . (4.20) 

Proof We set f{t,t') = 'D{c{t),c{t')) for t,t' G (-2e,2e). Using that D is symmetric 
and that D(x,x) = 1, we know hat 
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Thus the hnear term in a Taylor expansion vanishes, 

/(to, t) = ^ + \ 9{to) {t - tof + 0{\t- tof) , 

where we set 

As the function g is smooth, we can again expand it in a Taylor series, 

g{to)=g{0) + 0{to). 

We thus obtain 

fito,t) = 1 + ^ 5(0) (t - tof + O{\to\ \t - tol') +0[\t- tol') . 
Setting oi = 2(^(0) and using that \t\, |to| < 2e, the result follows. □ 

Lemma 4.22. Under the assumptions of Theorem \4 -151 the statement ()4.18p holds. 

Proof. Assume conversely that for x E supp/9 there is a point y £ suppp n /C(x). We 
choose the chain xq, xi = x, X2 as in Lemma [4. 211 We use the notation of Corollary l4.6l 
in case = 3, setting x^ = y. Choosing the vector u £ as u = (1, —2, 1, 0), we can 
apply Lemma 14.211 to obtain 

{u, Lu)c4 =6-4 V{xo, xi) + 2 a3(xo, X2) - 4 D(xi , xa) = O(e^) . 

Furthermore, using (j4.1ip . we know that 

Vixo,y) = be + Oie^) 

with 6 7^ 0. Thus, choosing u = {a, —2a, a, (3) with a, /3 G M, it is 

For sufficiently small e, the matrix in this equation has a negative determinant, in 
contradiction to Corollary 14.61 □ 

This completes the proof of Theorem 14.151 

In order to finish the proof of Theorem I4.17( we first remark that combining the 
symmetry of 2) with the assumption that T> is locally translation symmetric at x with 
respect to c, we know that T){c{t),c{t')) = f{\t — t'\). A Taylor expansion of / yields 
the following simplification and generalization of Lemma 14.211 

K 

D{c{t),c{t')) = l + Y^ai{t- t'f + 0({t- t')'('^+'^) , (4.21) 

i=l 

where the real coefficients Oj only depend on c. 

Lemma 4.23. Under the assumptions of Theorem \4-.n , the statement (I4.18P holds. 

Proof. Let us first verify that in the real analytic case, there is a p such that (I4.13P 
holds. Namely, assuming the contrary, all the t-derivatives of the function D(c(t),y) 
vanish. As the function T){c{t),y) is real analytic in a neighborhood of t = (as the 
composition of analytic functions is analytic), it follows that this function is locally 
constant. This contradicts the fact that 2)(c(t),?/) changes sign at t = 0. 
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Assume conversely that for x G supp p there is a point y G supppn/C(x). We choose 
the chain xo,xi = x,X2, . . . as in (j4.19p with k = p + 1. We use the notation of 
Coroharv 14.61 in case N = k. Then the Gram matrix L becomes 

/ 1 fie) ••• f{k8)\ 

/(e) 1 
i= (/(eN-j|)),,,=o,...,fc= : 

\f{ke) ••• 1 

Using the expansion (|4.21l) for K = k — 1, we obtain 

L = E + aie^{\i-jf)+a2e^{\i-j\^) 

+ ... + afe_ie2{fc-i) (I • _ j\2ik-i)^ ^ 0(^2fc 



(4.22) 



where E denotes the matrix where all the matrix entries (also the off-diagonal entries) 
are equal to one, and — jl"^) is the matrix whose element has the value |i — jl"^- 
Let us construct a vector v G C*'^^ such that the expectation value {v, Lv) is 0{e'^^). 
To this end, we take for v = {vi)\^Q G C'^^"'^ a non-trivial solution of the k linear 
equations 

k k k k 

Yvi = 0, Y,iv, = 0, Yi\, = 0, ...,J2i^'^Vi = 0. (4.23) 

i=Q i=0 i=0 i=Q 

Then {v, Ev) = and for alU G {1, . . . , /c - 1} 



J 3 

V 
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Each summand involves a power of i and a power of j, where always one of these 
powers is smaller than k. Thus all summands vanish according to (j4.23p . The solution 
V can always be normalized by t>o = 1, because setting to zero, the system of 
equations (I4.23P can be rewritten with the square Vandermonde matrix which has a 
trivial kernel. In view of the expansion (j4.22p . we conclude that {v,Lv) = 0{e'^^). 

We next consider the setting of Corollary 14.61 in case = A: + 1 and x^+i = y. 
Using (I4.13P together with the fact that the points y and xq are timelike separated, 
we find that 

C{xo,y) = beP + 0{eP+') (4.24) 
for 6 7^ 0. We choose the vector u G C'^"'"^ as n = {avQ, . . . , av^, (3) with a, /3 G M, 

where we combined (j4.24p with our normalization vq = 1, and used that y is not 
timelike separated from xi, . . . , x^. For sufficiently small e, the matrix in this equation 
has a negative determinant, in contradiction to Corollarv 14. 6i □ 



This completes the proof of Theorem 14.171 
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Figure 4. Numerical minima for the weighted counting measure on the circle. 



5. The Variational Principles on the Circle 

As a simple starting point for a more detailed analysis, we now consider the varia- 
tional principles on the circle (see Example (b) on page [5]). We first discuss numerical 
results, which again show the "critical behavior" discussed in Section [3] for the vari- 
ational principle on S'^. Applying the previous structural results, we will prove this 
critical behavior and show under generic assumptions that the minimizing measure is 
supported at a finite number of points. Moreover, we will give many minimizers in 
closed form. 

The numerical solution methods and results are similar as on S"^, as we now describe. 
We again consider the weighted counting measure (13. 3p . As the starting configuration 
we choose in analogy of the Tammes distribution on S"^ a uniform distribution of m 
points on the circle, 

{x, = e'^^-'^''-,k = 1, . . . , m} , = — , (5.1) 



X„ 



ni 



with uniform weights pk = 1/m. Minimizing as in Section [3] with a simulated annealing 
algorithm, we obtain the result shown in Figure HI The numerical results indicate that 
the minimizing measure is supported at a finite number of points mo. This number 
can be stated explicitly by 

mo = min <j n G N : n > — > , (5-2) 



where "dmax, as given by ()2.7p . denotes the opening angle of the lightcone. The number 
mo increases with r, with discontinuous "jumps" at the values 



1 - COs(t?m) ' 



(5.3) 



Besides the discrete nature of the minimizers, the numerical results reveal that at r = 
Tc = \/2 (corresponding to t?max = f )> the structure of the minimizers changes com- 
pletely. Just as in Section [3l this effect can be understood as a phase transition. More 
precisely, if r < Tc, every minimizer is generically timelike. If we further decrease r 
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(i.e.. for every fixed r < t^), we even found a large number of minimizing measures, 
supported at different numbers of points with strikingly different positions. However, 
if r > \/2, the minimizer is unique (up to rotations on S^), is supported at uiq points, 
and is not generically timelike. 

In the remainder of this section, we make this picture rigorous. First, the opera- 
tor can be diagonalized explicitly by plane waves 4>n{x) = e*"'^"' (where n G Z, 
and -dx is the angle). This gives rise to the decomposition 

2 

n=l 

where 

1^0 = [[ Dix, y) di^ix) dfi{y) = At^ - . (5.4) 

and similarly ui = 2t^ and 1^2 = ^t^- In the case r < 2 all eigenvalues uq, vi and 1^2 
are non-negative, and we can apply Proposition 14.111 to obtain 

^min ^ '^O • 

For sufficiently small r, the uniform distribution of points on the circle (15. ip gives a 
family of generically timelike minimizers. 

Lemma 5.1. If m > 3 and r is so small that C{x,y) = T){x,y) for all x,y £ Xm, 

then p = ^ X^I^i is a generically timelike minimizer. Every other minimizer is 
also generically timelike. 

Proof. A straightforward calculation using the identities 

m—1 m—1 ^ 

^ e^'^'^- =0 and ^ 
fc=0 fc=o 

yields for any x £ S^, 

^ m—1 

d(x) = — V (2 + 2(X, Xk) - + T^ix, Xkf) 

m ^-^ 

k=0 

^ - 2 /r,_ _ 2 , 1^ Jl 



-2t^ 2m - mr^ H = i/q • 

m \ 2 / 

In particular, one sees that S[p] = uq. 

The assumption jC{x,y) = D(x,y) for all x,y £ Xm can only be satisfied if r < 2. 
Thus in view of ()5.4p . the operator is positive semi-definite. We finally apply 
Proposition 14.111 □ 

Applying this lemma in the case m = 4 gives the following result. 

Corollary 5.2. If t < Tc, every minimizer is generically timelike. 

More general classes of generically timelike minimizers can be constructed explicitly 
with the help of Corollarv l4.12[ In particular, one can find minimizing measures which 
are not discrete. For the details we refer to the analogous measure on S"^ given in 
Example [ 
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Figure 5. A minimizer for r = 4. 



Having explored the case r < Tc, we proceed with the case r > Tc. As already 
stated, the closed lightcones are given by 

Jix) = |y : (x, y) > 1 - = cos(t?max)| U {-x} . 

Therefore if r > \/2 = (or equivalently i^max < f )) the condition of antipodal points 
(see Proposition 14.13]) is satisfied. Thus there are no generically timelike minimizers. 
As the condition (14. lip is obvious, we can apply Theorem 14.151 (A) and conclude that 

if r > Tc, every minimizing measure is discrete . (5-5) 

Using results and methods from Section [4. 3^ we we will be able to explicitly construct 
all minimizers under the additional technical assumption that 



r > Td := y3 + \/lO. 
We first introduce a descriptive notation: 

Definition 5.3. A chain of length k is a sequence xi, . . . , € of pairwise distinct 
points such that (xj, Xj+i) = cos(T?max) for all i = 1, . . . , k — 1. 

Theorem 5.4. If t > t^, the support of every minimizer p is a chain {xi, . . . ,Xmo} 
(with mo as given by (15. 2p ). The minimal action is 

£(0)(£(0)+£(7)) 



(mo-2)(£(0)+£(7)) + 2£(0) ' 



(5.6) 



where 7 = arccos((xi, Xmo)) ^ (0,'i?max]- The minimizing measure is unique up to 
rotations on S^. 

An example for the support of the minimizing measure is shown in Figure El Up to 
rotations, the points of the chain can be written as 

Xk = e^('=-^)'^-- , k = l,...,mo. (5.7) 

In the special cases t = Tm, the minimizer is the measure with equal weights supported 
on the uniform distribution X„i- In the general case, the weights will not all be the 
same, as will be specified below. 

For the proof of Theorem 15.41 we proceed in several steps. 

Lemma 5.5. If t > ^/6, the minimal action is attained for a measure supported on a 
chain xi, . . . ,Xk- In the case k = mo, every minimizing measure is a chain. 
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Proof. Let /? be a minimizing measm'e. We first note that every chain K in the support 
of p must have finite length, because otherwise i^max/^r would have to be irrational. 
As a consequence, K would be a dense set of S^, in contradiction to the discreteness 
of p (see (jS.Sp ). Let us assume that the support of p is not a chain. 

We let K C suppp be a chain, which is maximal in the sense that it cannot be 
extended. Set L = suppp \ K. We consider variations of p where we rotate K hy a 
small angle i), leaving the weights on K as well as p\l unchanged. The fact that K 
cannot be extended implies that that these variations are smooth in at = 0. The 
minimality of p implies that 

6S = and 6^S = ^ 2 p{x) p{y) 6^ C^x , y) > . (5.8) 

xGK,y&L 

On the other hand, differentiating (j2.6p . one finds that the function D restricted 
to [0, '& 

maxj IS concave, 

V"{^) = -4r2(cos(??) + cos(2t9)) < {if t > Vg) . (5.9) 

Comparing with (j5.8p . we conclude that C{x,y) vanishes for all x £ K and y £ L. In 
the case that i^K = niQ, this implies that L = 0, a contradiction. In the remaining 
case < mo, we can subdivide the circle into two disjoint arcs Ak and Al such 
that K C Ak and L C Al. The opening angle of Ak can be chosen larger than "i^max 
times the length of K, giving an a-priori upper bound on the length of K. 

By further rotating K, we can arrange that the chain K can be extended by a point 
in L, without changing the action. If the extended chain equals the support of p, the 
proof is finished. Otherwise, we repeat the above argument with K replaced by its 
extension. In view of our a-priori bound on the length of K, this process ends after a 
finite number of steps. □ 



Lemma 5.6. Suppose that p is a minimizing measure supported on a chain. If t > 
Vo + VTO, the length of this chain is at most mo. 

Proof For all 7 G (0, 

max) 3-11 elementary calculation shows that 

£(7)' + /:(T?max-7)' > AO)'. (5.10) 

In the case r = there is nothing to prove. Thus we can assume that r / Tm^- 
For a chain xi, . . . ,Xk with k > mo, the Gram matrix corresponding to the points 
xi, Xmo+i, X2 has the form 

£(0) /:(T?max-7) \ 

A^max-7) m At) • (5.11) 

ml 

Using ()5.10p . its determinant is negative, in contradiction to Corollarv 14.61 □ 

From the last two lemmas we conclude that every minimizer p is supported on one 
chain of length at most mo- Parametrizing the points as in (j5.7p . the only contributions 
to the action come from C{xi,xi) and C{xi,Xmo)- Using Lagrange multipliers, the 
optimal weights pi = p{xi) are calculated to be 

i , _^ 

£(0)+£(7) ^^"£(0) 



Pi = Pmo = rrn\ , rf.^ Pi = 7v7^ for i = 2, . . . , mo - 1 , (5.12) 
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where we set 

£(0)(£(0)+£(7)) 
(mo-2)(£(0) + £(7)) + 2£(0) ' 

The corresponding action is computed to be S[p] = A, giving the formula in (15.6p . 
Using this expHcit value of the action, we obtain the fohowing 

Lemma 5.7. Suppose that p is a minimizing measure supported on a chain. Then the 
length of this chain is at least mo- 

Proof. For a chain of length n < mo, the only contribution to the action come 
from X/), / = 1, . . . , n. The corresponding optimal weights are computed by pi = 
1/n. The resulting action is 

" 1 1 
5 = -^C{xi,Xi) = -C{0) . 

This is easily verified to be strictly larger than the value of the action in (15. 6p . □ 

This completes the proof of Theorem 15.41 

We finally remark that if r lies in the interval (\/2, \/3 + \/lO) where Theorem 15.41 
does not apply, the numerics show that the minimizing p is again the measure sup- 
ported on the chain of length mo, with one exception: If r is in the interval (1.61988, rs) 

with Ts = ^2 + a chain of length mo + 1 = 6 gives a lower action than the chain 

of length 5. In this case, the Gram matrix ()5.1ip is indeed positive definite, so that 
the argument in Lemma 1 5 . 71 fails . 



6. The Variational Principles on the Sphere 

We now come to the analysis of the variational principles on the sphere (see Exam- 
ple (a) on page 4). Applying Theorem 14.151 (A) with the curve c chosen as the grand 
circle joining x and y, we immediately obtain that every minimizing measure p on S"^ 
is either generically timelike or supp p = 0. The numerics in Section [3] indicated that 
these two cases are separated by a "phase transition" at r = Tc = \/2. We will now 
prove that this phase transition really occurs. Moreover, we will develop methods for 
estimating the minimal action from above and below. Many of these methods apply 
just as well to the general setting introduced in Section [2] (see (I2.ip - (l2.5p ). 

6.1. Generically Timelike Minimizers. We first decompose 2) in spherical har- 
monics. A short calculation yields in analogy to (14. 7p the decomposition 

2 I 

a)(x,y)=i/o + 47r^j.z ^ l^™(x)lf^, 

1=1 m=—l 

where the eigenvalues are given by 

9 4 4 4 4 . , , 

^/o = 4T2--r^ i^i = -t\ ^2 = — • (6.1) 
3 o 15 

In particular, the operator is positive semi-definite if r < -v/3. 

If T < Tc, there is a large family of minimizers, as we now discuss. The simplest 
example is the octahedron: Denoting the unit vectors in by 61,62,63, we consider 
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the measure p supported at itej with equal weights g. Obviously, the condition (i) in 
Definition 14.71 is satisfied. Moreover, for any x E one calculates 

A{x) = \ 2T2(2 + 2(x,y)-r2 + r2(x,y)2) = 

?/6suppp 

= (12 - 6r2 + 2t\xI + xl+ xD) = . 

Thus Proposition 14.111 yields that /? is a generically timelike minimizer. Moreover, 
from Proposition 14. 1 1] we conclude that every minimizer is generically timelike. If con- 
versely T > Tc, the condition of antipodal points is fulfilled, and thus Proposition 14.131 
shows that no generically timelike minimizers exist. We have thus proved the following 
result. 

Corollary 6.1. If t < Tc, every minimizing measure p on 5*^ is generically timelike, 
and the minimal action is equal to vq as given by (j6.ip . // conversely r > Tc, every 
minimizing measure p is not generically timelike and supp p = 0. 

Using Corollary 14.121 one can also construct minimizers which are not discrete, as 
is illustrated by the following example. 

Example 6.2. We introduce the function / G L^(S'^) by 



''I if G [0,arccos(0.8)], 

^ if ?? G [arccos(0.4), arccos(0.2)] 

^ if ?? G [arccos(— 0.5), arccos(— 0.7)], 

otherwise. 



Then if r < 1.00157, a straightforward calculation shows that / has the properties (a) 
and (b) of Corollary 14.121 Thus the measure dp = fdp is a minimizing generically 

timelike measure with supp p / 0. 

6.2. Estimates of the Action. As not even the solution of the Tammes problem is 
explicitly known, we cannot expect to find explicit minimizers for general r. Therefore, 
we need good estimates of the action from above and below. We now explain different 
methods for getting estimates, which are all compiled in Figure [6l 

Estimates from above can be obtained simply by computing the action for suitable 
test measures. For example, the action of the normalized volume measure is 

S[p] = ^ sini? V{^) = 4 - ^ > . 

As one sees in Figure [6l this estimate is good if r is close to one. Another example is 
to take the measure supported at the Tammes distribution for K points, with equal 
weights. We denote the corresponding action by . We then obtain the estimate 

^min ^ St '■= min57^ . 

K 

One method is to compute St numerically using the tables in [14]. This gives quite 
good results (see Figure [6]), with the obvious disadvantage that the estimate is not 
given in closed form. Moreover, the Tammes distribution is useful for analyzing the 
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Figure 6. Estimates of the action on S^: Upper bounds obtained from 
the volume measure S[fj] and from the Tammes distribution St, lower 
bounds by uq and by the heat kernel estimate Sk- 



asymptotics for large r. To this end, for every Tammes-distribution Xk we intro- 
duce tk as the minimal value of r for which all distinct points in Xk are spacelike 
separated. In analogy to (|5.3p . the value of tk is given by 



TK 



1 — cos{'dK) 



where i^k now denotes the minimal angle between the points of the Tammes distribu- 
tion, 

T^K = max min arccos((xj, Xj)) . 

Using an estimate by W. Habicht and B.L. van der Waerden for the solution i^k 
(see [13l page 6]), we obtain 



r^V'^^V%4>4^ (6.2) 



for some constant C > 0. For given r > 1 we choose K gN such that tk-i < t < tk- 
Then 

c cK-i_ 8r2 84 K //8^\^/2 C ^ 



In the limit r — )• cxd, we know that K ^ oo, and thus 

4^/3 



limsup5min > 



Constructing a lower bound is more difficult. From (j6.ip it is obvious that the 
operator T>^ is positive semi-definite if r < \/3- Thus we can apply Proposition 14.111 
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to obtain 

5min > z^o if r < \/3 . 
If r < -v/2, this lower bound is even equal to 5min according to Corollary 16. 11 As shown 
in Figure [6l the estimate is no longer optimal if r > \/2- 

Another method to obtain lower bounds is based on the following observation: 

Proposition 6.3. Assume that is an integral operator on with integral ker- 
nel K G C^{S^ X S'^jM) with the following properties: 

(a) K{x, y) < C{x, y) for all x,y £ S"^. 

(b) The operator is positive semi- definite. 
Then the minimal action satisfies the estimate 

K{x,y) d^i{x)d^i{y) . (6.3) 

Proof. For any p S assumption (a) gives rise to the estimate 

S[p]=ll C{x,y)dp{x)dp{y) > K{x,y)dp{x)dp{y) . 

Next, using property (b), we can apply Proposition 14.31 to conclude that the volume 
measure is a minimizer of the variational principle corresponding to K, i.e. 



K{x, y) dp{x) dp{y) > K{x, y) dp{x) dp{y) . 

S^xS^ JJs^xS^ 

Combining these inequalities gives the result. □ 
In order to construct a suitable kernel, we first consider the heat kernel ht on S'^, 

oo I 

ht{x,y) = (e*^s2) (x,y) =4^^ e-* J] Yr{x)Yr¥) . 

1=0 m=-l 

The heat kernel has the advantage that condition (b) is satisfied, but condition (a) is 
violated. This leads us to choosing K as the difference of two heat kernels, 

K{x, y) = X {ht, {x, y) - 6ht, {x, y)) . (6.4) 

For given ti < t2, we choose 6 and A such that K{x,x) = 1 and i^('!9max) = 0, i.e. 

= — tt; r < 1 and A = - — — — -— > . 

By direct inspection one verifies that condition (a) is satisfied (see Figure[7]for a typical 
example). The eigenvalues of the operator are computed to be 

A(e-*i'a+i)-5e-*2'('+^)), 

showing that the operator Kf^ is indeed positive semi-definite. Thus we can apply 
Proposition 16.31 Using that 



ht{x, y) dfi{x) dp{y) = Att Y^{x)Yo^{y) dp{x) dfi{y) = 1 

52x52 JJs^xS^ 

we obtain the heat kernel estimate 

^min > 5"/^ = A (1 - (5) . 
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Figure 7. The Lagrangian L and the function K in the heat kernel 
estimate for r = 2. 

In this estimate, we are stih free to choose the parameters t\ and t^- By adjusting 
these parameters, one gets the lower bound shown in Figure O Thus the heat kernel 
estimate differs from the minimal action only by an error of about 20%, and describes 
the qualitative dependence on r quite well. But of course, it does not take into account 
the discreteness of the minimizers. 

7. The Variational Principles on the Flag Manifold 3"^'^(C-^) 

We finally make a few comments on the variational principles on the flag mani- 
fold ji'2(C-'') (see Example (c) on page Ej) . We first apply our main Theorems 14.151 
and 14.171 to obtain the following general result. 

Theorem 7.1. Every minimizer p on 3""*^'^ is either generically timelike or supp p = 0. 

Proof. As a homogeneous space, the flag manifold 3^^''^{C^ ) has a real analytic structure 
(see [111 Chapter II, §4]). Then the function D is obviously real analytic. Moreover, 
it is symmetric and constant on the diagonal. In order to apply Theorem 14.151 for 
given y S /C(x) we must find a curve c joining x and y which satisfies ()4.1ip . Al- 
ternatively, in order to apply Theorem 14.171 our task is to construct a curve c(i) 
with c(0) = X and c(l) = y which is analytic in a neighborhood of t = 0, such that the 
function D{c{t),y) changes sign at t = 0. Since in the case r = 1, the sets }C{x) are all 
empty, we may assume that r > 1. 

We denote the range of a; by I C C-^ and the orthogonal projection to / by ttj. 
Choosing an orthonormal basis (61,62) of I, the matrix x\i can be represented with 
Pauli matrices by 
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where "Q denotes the angle between u and v. The operator y has rank two if and 
only \i K > \p\. A short calculation shows that in this case, 2) only has transverse 
zeros. Thus we can choose a direction c(0) where the condition (j4.1ip is satisfied. 
Choosing a smooth curve starting in this direction which joins x and y, we can apply 
Theorem 14.151 (A) to conclude the proof in this case. 

It remains to consider the situation when y has rank at most one. This leads us 
to several cases. We begin with the case when y\i vanishes. In this case, we may 
restrict attention to the four-dimensional subspace U = Imx © Imy. In a suitable ba- 
sis (ei, . . . , 64) of this subspace, the operators x and y have the matrix representations 



where again u,v £ S^. A unitary transformation of the basis vectors ei and 62 describes 
a rotation of the vector -u in R^. By a suitable transformation of this type, we can 
arrange that the angle between u and v equals ??max (see (12. 7p ). We now define the 
curve c : [0, vr] 3"^'^ by 



where w : [0, vr] — 5^ is the geodesic on 5^ with w{0) = u and 'w{7r) = v. The curve c 
is a real analytic function with c(0) = x and c{tt) = y, which is obviously translation 
symmetric. Furthermore, one computes 



where T)g2 is the corresponding function on the unit sphere ()2.6p . As 1)52(1?) changes 
sign at t?max5 the function D(c(t),y) changes sign at t = 0. Thus Theorem 14.171 (A) 
applies, completing the proof in the case y\j = 0. 

We next consider the case that y has rank one. We choose the basis (61,62) of / 
such that y is diagonal. 



An elementary consideration shows that we can extend the basis of / to an orthonormal 
system (61,62,63) such that on the subspace J := ({61,62,63}), the operator y := 
TTjyirj has the form 





(7.1) 



V{c{t),y) = sm{t)^ 'Ds2{w{t),v) 




with a 7^ 0. 




with a ^ and ac 7^ |6p . 



(7.2) 




(7.3) 



Setting y{t) = U{t)yU{t) the matrix y becomes 




where p and k are new parameters with 

K > \p\ and 



p + KVs = a^O 



(7.4) 
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and V G S'^ is again a unit vector. The function T> is now computed by 
^i^,y{t)) = ^Tr(x|7y(t))'-2det {x\iy{t)) 

= i (Tr{x\j yit)f - 4 (t2 - 1) (k^ - p^) sin^ t) . (7.5) 

In order to simplify the trace, we transform the phase of 63. This changes the phase 
of b in ()7.2p . thus describing a rotation of the vector v in the (1, 2)-plane. This makes 
it possible to arrange that the vectors (ni,n2) and (^1,^2) are orthogonal in M^. We 
thus obtain 

TT{x\iy{t)) = (1 + TU3){p + KV3) + (1 - TU3){p - KV3) sin^t . 

We now have two subcases: 

(1) V3 7^ ±1: We vary the vectors u and v as functions of t such that the above 
orthogonality relations remain valid and 

U3 = cos{'& + at) , V3 = cos{ip + f3t) 

with free "velocities" a and /?. Since C{x, y) = at t = 0, we know that 

1 v/r2 - 1 

cos'& = , sin'&=- /O. (7.6) 

T r 

A Taylor expansion yields 

Tr(x|/y(t)) = —t ar {p + KV3) sint? (7.7) 

+ — (^(4 + a^) p + (-4 + a^)K cos (/? + 2af3KT sin 'd sin v^j + 0{t^) . (7.8) 

As the factor (p + ^^3) is non-zero in view of (j7.4p . the linear term (j7.7p does 
not vanish whenever a 7^ 0. By suitably adjusting a, we can arrange that the 
square of this linear term compensates the last term in (j7.5p (which is also 
non-zero in view of our assumption k > p). Next, we know from (j7.6p and 
our assumptions that the term ~ a(3 in (j7.8p is non-zero. Thus by a suitable 
choice of f3, we can give the quadratic term (|7.8p any value we want. Taking 
the square, in (j7.5p we get a contribution ~ t^. Thus the function D changes 
sign. Transforming to a suitable basis where y is a fixed matrix, we obtain a 
curve x{t) which is locally translation symmetric. Extending this curve to a 
smooth curve c which joins the point y, we can apply Theorem 14.171 (A). 

(2) U3 = ±1: We know that the matrix y is diagonal, 

Now we keep v fixed, while we choose the curve u{t) to be a great circle which 
is inclined to the (l,3)-plane by an angle 7 / 0, i.e. 



U3 = cos(t? + at) cos 7 . 



Repeating the above calculation leading to (j7.7p and (|7.8p . one sees that we 
again get a non-zero contribution to D of the order ~ t"^. Thus D again changes 
sign, making it possible to apply Theorem 14.171 (A). 
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It remains to consider the case when y vanishes but y\i ^ 0. A short consider- 
ation shows that y\i cannot have rank two. Thus we can choose the orthonormal 
basis (61,62) of / such that yei ^ and 2/62 = 0. By suitably extending this or- 
thonormal system by 63 and 64, we can arrange that the operator y is invariant on the 
subspace ({61,62,63,64}) and has the matrix representation 



2/|({ei,e2,e3,e4}) 



/O 

a 

\0 





b 



If 6 7^ 0, we can again work with the curve ()7.ip . If on the other hand b 
operator y is invariant on ({61,62,63}) and has the canonical form 



0, the 



y\{{ei,e2,ea}) 








Transforming y by the unitary matrix 

V{t) 



cos T 

- sinr 



smr 
cos T 





2 



we can arrange that y is again of the form (17. 2p . but now with coefficients depending 
on r. Setting t = t'^, we can again use the construction after (17. 2p . This completes 
the proof. □ 

For sufficiently large r, we can rule out one of the cases in Theorem 17.11 showing 
that the minimizing measures do have a singular support. 

Theorem 7.2. There are no generically timelike minimizers if 



3/ + 273(F 



1 



(2 + /) 

The method of proof is to apply Proposition 14. 13] (I) . In the next two lemmas we verify 
the necessary assumptions and compute i'q. 

Lemma 7.3. The operator has rank at most 3/^. 

Proof. We extend the method used in the proof of [Sj Lemma 1.10]. A point x G 9" is 
a Hermitian / x /-matrix of rank two, with non-trivial eigenvalues 1 + r and 1 — r. 
Thus, we can represent x in in bra/ket notation as 

X = \u{x)){u{x)\ — \v{x)){v{x)\ , 

where u{x) and v{x) are the eigenvectors of x, normalized such that 

{u{x)\u{x)) = T + 1 and {v{x)\v{x)) = r — 1 . 

A short calculation shows that the non-trivial eigenvalues of the matrix product xy 
coincide with the eigenvalues of the 2 x 2-matrix product 

'{u{x)\u{y)) -{u{x)\v{y))\ f{u{y)\u{x)) 
^{v{x)\u{y)) -{v{x)\v{y)) ) \{v{y)\u{x)) 

Using (|2.8p . we can thus write the function T> as 

, 2" 



A 



-{u{y)\v{x)) 
-{v{y)\v{x)) 



T>{x, y) = Tr 



A 



xy 
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This makes it possible to recover T){x,y) as the "expectation value" 



V 1^ V* (gi V (gi V* 




of a suitable matrix B, whose 3x3 block entries are of the form 

Bij = bij + 5i^2Sj,2 (cipi + C2P2 + C3P3) with bij,Ci G 
and the operators pi permute the factors of the tensor product, 



pi(u<^u 

P2 {U (g) U* 
P3 {u (g) U* 

Hence introducing the operator 



' V (g V ) = V (g> V (giU(g>U 
I V V*) = U <g V* V <g u* 

'V(E>v* ) = V u* (E> u g) V* 



U(g) u (g u(g> u 
u®u*®v®v* 




V ig V (XI V ig V 

we find that D^j = K*BK. This gives the claim. □ 

In view of this lemma, we may decompose D in the form ()4.7p . 

Lemma 7.4. The eigenvalue vq in the decomposition (j4.7p is given by 

2(3/ + 6/T^-(2 + /)r4-6) 
-0 = j^j2^^ . (7.10) 

Proof. Now it is now most convenient to represent the elements in 3~ as 

{1 + t)\u){u\ + {1-t)\v){v\ , (7.11) 

where the vectors n, u G C'^ are orthonormal. Then the normalized volume measure ^ 
on 3" can be written as 



dp 



1 



6{Re {u,v)) 6{lm {u,v)) d{\\uf - l) S{\\vf - l) dudv, 



vol(3") 

where du and dv denotes the Lebesgue measure on C-^. The total volume is computed 
to be 



vol(3") 



5{Re {u,v))5{lm {u,v))5{\\u\\ 



1 



CfxCf 
^V0l(52^~l)v0l(S2/-3) 



f 11^ — l)du dv 



To simplify the calculations, we fix x and choose an eigenvector basis of x. Then x = 
diag((l + r), (1 — r), 0, . . . , 0), whereas y is again represented in the form ()7.1ip . Then 
the eigenvalues of the product xy only depend on the components ui,U2,vi,V2- More 
precisely, using ()2.8p . we obtain 



V{x,y) =i[(l + r)2|niP + (l-r2)(|z;i|2 



F2I 



[l-rf\v2\'f 



+ 2(1 - T^)! (1 + r)uiul + (1 - T)viv^f =: f{u, v). 
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Our task is to compute the integral z/q = f{u, v) dfi. In the case / > 4, one uses the 
symmetries to reduce to a lower-dimensional integral, 

foo foo roo roo f f /"OO 

uo = c dui / du2 / du'i I dvi / dv2 / dv^ / dv4 
Jo Jo Jo Jo Jc Jc Jo 

X 6{\\uf - 1) 6{\\vf - 1) 6 (Re ((n, v))) 5 (Im {{u, v))))fiu, v) uxU2v't\xv't\ 
where c is the constant 

1 



vol(52^-^) vol(S2-^-^) (27r)3 . 



vol(J-) 

Now carrying out all integrals gives the claim. The proof in the case / = 3 is similar. 

□ 

The remaining question is whether generically timelike minimizers exist for small r. 
In the special case r = 1, the operator 2)^ = is positive semi-definite (see jH 
Lemma 1.10]), so that Proposition 14.31 or similarly Proposition 14.111 vields that the 
standard volume measure is a generically timelike minimizer. However, if r > 1, 
these propositions can no longer be used, because the operator fails to be positive 
semi-definite: 

Lemma 7.5. // / > 3 and r > 1, the operator has negative eigenvalues. 

Proof. Since supp/i = 3", it suffices to find two points xi,X2 G 3" such that the corre- 
sponding Gram matrix T>(xi,Xj) is not positive semi-definite. For given e S (0, 1) we 
choose the four vectors 

ui = ei , vi = 62 and n2 = ei , V2 = Ve €2 + \/l - e 63 

(where e, are the standard basis vectors of C-^) . Taking the representation (j7.1ip , we 
obtain two points xi,X2 G 3". The corresponding Gram matrix is computed to be 

( H-<r-l? + {r + irf\ 

[U-e{r - If + (r + Iff 8r^ J' 

The determinant of this matrix is negative for small e > 0. □ 

In this situation. Proposition 14. 14] still gives some information on the possible support 
of generically timelike minimizers. However, it remains an open problem whether and 
under which conditions generically timelike minimizers exist. 
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